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O ■ Abstract 

; 

We study the dynamics of a quantum superconducting circuit which is the 

c/3 ' analogue of an atom in a high-Q cavity. The circuit consists of a Josephson 
charge qubit coupled to a superconducting resonator. The charge qubit can 
be treated as a two level quantum system whose energy separation is split by 

1 , the Josephson energy Ej. The superconducting resonator in our proposal is 
^ , the analogue of a photon box and is described by a quantum harmonic oscilla- 
O \ tor with characteristic frequency u^r- The coupling between the charge qubit 
' and the resonator is realized by a coupling capacitance Cc. We have studied 

^ ■ the eigenstates as well as the dynamics of the quantum circuit. Interesting 

^ , phenomena occur when the Josephson energy equals the oscillator frequency, 

*^ \ Ej = hujj.. Then the quantum circuit is described by entangled states. We 

I have deduced the time evolution of these states in the limit of weak coupling 

OO ■ between the charge qubit and the resonator. We found Rabi oscillations of 

^ ■ the excited charge qubit eigenstate. This effect is explained by the sponta- 

, neous emission and re-absorption of a single photon in the superconducting 

"ti ! resonator. 

■ I. INTRODUCTION 

c : 

Cj . Recently, a substantial interest in the theory of quantum information and computing0 

^ ! has revived the physical research on quantum systems. The elementary unit of quantum 

\ information is a two-state system, usually referred to as a quantum bit (qubit). Basic 
' operations are realized by preparation and manipulation of, as well as a measurement on, 
entangled states in systems which consist of several coupled qubits. However, the fabrication 
of physical systems which would enable the actual implementation of quantum algorithms is 
far from being realized in the near future and a substantial amount of fundamental research 
is still needed. 

During the past five years, great progress has been made in the manipulation of entangled 
states in systems consisting of up to four qubits based on ion trapJ^ and atoms in a high-Q 
cavityi. These two experiments demonstrate clearly and unambiguously the possibility to 
coherently control the entangled states of a limited number of qubits, as well as to perform 
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a quantum measurement on them. In spite of this success, it seems quite difficult to reahze 
circuits consisting of the large number of ion traps or atoms in a cavity necessary for quantum 
computation. 

It has been suggested that small solid state devices fabricated using nanolithography 
technologies are promising for quantum circuit integration. However, the coherent ma- 
nipulation of entangled states as well as the realization of quantum measurements remain 
fundamental issues to be investigated. One of the main challenges is to gain control over all 
possible sources of decoherence. At present the best candidates for the implementation of 
quantum gates based on solid state devices are circuits using small Josephson junctions. In 
the superconducting state, such circuits contain less intrinsic sources of decoherence. Indeed 
it has been experimentally demonstrated that a single Cooper pair box is a macroscopic two 
level system which can be coherently controlledoQ. At about the same time, theoretical 
works have proposed the use of the Cooper pair box as a qubit (the so-called Josephson 
charge qubit) in the context of quantum computers. In particular, systems consisting of 
several charge qubits with controlled couplings have been discussed, the quantum measure- 
ment problem has been addressed and the decoherence time has been More 
recently, qubits based on superconductingloops containing small Joseohson junctions have 
been proposed (Josephson phase qubits)li0 and are currently studied0'0. But up to now, 
the existence of entangled states, which are at the heart of quantum information processing, 
has been demonstrated neither for charge qubits nor for phase qubits. 

In this article we propose to study one of the simplest Josephson circuits in which en- 
tangled states can be realized. It consists of a charge qubit coupled to a superconducting 
resonator and can be described theoretically by a two level system coupled to a quantum 
harmonic oscillator. After a description of this quantum circuit in the next section, the 
Hamiltonian describing it will be derived in Sec. |in[ In Sec. fV|, the time evolution of the 
eigenstates is obtained and we demonstrate the existance of entanglement. In the last Sec- 
tion, we discuss the dynamics of the quantum circuit for typical experimental values of the 
system parameters. 



II. QUANTUM CIRCUIT 

The Josephson circuit we study hereafter is depicted in Fig. |1]. It consists of three different 
elementary circuits: a Cooper pair box, an LC-resonator and a "coupling" capacitor. 

For small enough junction capacitance Cj, gate capacitance Cp, and coupling capacitance 
Cc, the charging energy of the box is large compared to thermal fluctuations and the excess 
charge of the box is quantized. On the other hand, we assume the charging energy to be 
smaller than the superconducting gap A, such that no quasiparticles are present in the box. 
Thus the excess charge is entirely due to the presence of Cooper pairs and charge quantization 
occurs in units of 2e. The gate voltage Vg is used as an external control parameter. When 
the gate charge Ng = —CgVg/e is equal to unity, the Cooper pair box can be viewed as a 
macroscopic two-level quantum system whose energy separation is split by the Josephson 
energy Ej. The two eigenstates |— ) and |+) correspond to a coherent superposition of the 
two different charge states of the boxH. When Ng ^ 1, the Cooper pair box will be referred 
to as a Josephson charge qubit. 
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The LC-resonator system can be described by a quantum harmonic oscillator whose 
characteristic frequency is given by uo^. This system is the analogue of a high-Q cavity. 

The capacitance Cc plays a crucial role in our proposed circuit since it couples the charge 
qubit and the resonator to each other. These two circuits are no longer independent and 
the system must be considered in its totality. Thus the proposed quantum circuit of Fig. |1| 
realizes the simple situation in which a two level system is coupled to a harmonic oscillator. 
In spite of its simplicity, such a system describes a great variety of interesting situationsEHll. 



III. HAMILTONIAN 

The circuit depicted in Fig. |l| can be characterized mechanically by two generalized 
coordinates, (pj and (pr- These coordinates are associated with the voltage drop 5Vj over 
the junction and 5Vt over the resonator, respectively, according to the Josephson relation 
(pi = 2eSVit/h {i = j,r). We seek the Lagrangian C{(pj,(pr,(pj,(j)r) = T — V describing the 
dynamics of these variables. The potential energy V is a function of the coordinates only, 

V{(pj, (pr) = -Ej COS (pj + -^(pl, (1) 

where = (l/Lj,)(^/2e)^ is the energy associated with the inductance Lj. of the resonator. 
The kinetic energy T is quadratic in the velocities (pj and (pr- It is just the free electrostatic 
energy stored in the capacitators present in the circuit. This free energy can be written as 
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T = - 
2 



C^j{h(Pj/2ey + Cj^r{h(pr/2ey + 2C^,{h/2ey(P,(Pr . (2) 



Here we introduced the capacitances Cec = Cc + Cg, Csj = Cj + Csc, Csr = Cr + Cy,c- 
Note that the Lagrangian contains an interaction between the resonator and the junction: 
i^int = CEc(^/2e)^0j0r- The effective coupling between these two parts of the circuit is 
determined by the sum of the gate capacitance and the coupling capacitance. We also note 
that the equations of motion d{dC/ dcpi) / dt + dC/dcpi = express current conservation in 
the circuit. 

Through the Josephson junction. Cooper pairs can tunnel from or onto the island. The 
number of excess Cooper pairs on the island, n, depends on the gate voltage. Charge 
neutrality leads us to relation 

2ne = Cj^j{h(pj/2e) + C^c{'h(pr/2e) + NgC. (3) 

It is always possible to add a term to the Lagrangian which is a total time derivative. Let 
us add the term h(pjNg/2. As a result, the momenta conjugate to (pj and (pr are 

Pj = dC/d(pj = {h/2e)[Cj:j{h(pj/2e) + Csc{Hr/2e) + Ngc] = hn, 
Pr = dC/d(pr = C^r{h/2ef4)r + C^ci^i / 2ef 4) j . 

Note in particular that the momentum pj is proportional to n. 

The Hamiltonian is obtained with the help of the Legendre transform H = pj(pj+pr(pr—£^- 
We find H = Hj + Hr + He, where 
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Hj = Ecji2n - N,)' - Ej cos <f>j, (4) 

Hr^Ec,r{^Pr/nf + Er(t>l/2, (5) 

H,^-Ec,c{2n-Ng){2pr/n). (6) 
The charging energies Ec,i {i = j,r,c) appearing here are given by Ec,i = e^/2Ci^eff, with 

Q,eff = Cj + (1/Csc + l/Cr)-\ Cr,eS = + (l/Csc + l/Q)"', and C^.eff = [C^c +\l/Cj + 

l/Cr)~^][{Cj + Cr)/2C-£c]- The Hamiltonian equations of motion, pi — —dH/d(f)i, — 
dH/ dpi lead us again to current conservation. 

In order to obtain the quantum mechanical Hamiltonian H, wc replace Pi, 0j by the 
corresponding operators, with \pk,4>m] = ip-/i)^k,m- In particular, in 0— representation we 
have pk — {h/i)d/d(t)k- Note also that [n, (f)j\ — —i and n — —id/d(j)j. Below we discuss the 
various contributions to H in some detail. 

Josephson junction. The commutation relation [tt-, 0j] = implies [n, e**^^] = e**^^ . Using 
the basis states \n), where n corresponds to the number of excess Cooper pairs on the island, 
we thus have e''^^\n) = \n+ 1). Similarly, e~^'^^\n) = |n — 1). Therefore we can write Hj as 

Hj^Ec,jj:{2n-N,f\n){n\ 

n 

E(l^ + i)H + l^-i)H)- (7) 

^ n 

If the gate- voltage is such that A^^ ^ 1, the states with n = and n = 1 are almost 
degenerate. At low temperatures, the Hamiltonian Hj involves only these two states, and 
thus can be presented as a matrix 

This matrix can be diagonahzed. The eigenvalues are 

E^ = EcAl + m,f] T \^{5E,Y + E], (9) 

where 5Ng = Ng — 1 and 6Eg — —AEcjSNg. The corresponding eigenstates are 

|-) = «|0)+/3|1) (10) 
|+)=/3|0)-a|l) (11) 



where ^ 1 - = [1 + 5Eg/^{5EgY + E]]/2. 

Resonator. Since the LC-circuit constitutes just a harmonic oscillator with a charac- 
teristic frequency ujr — ^Jl/ LrCr,eS-i the Hamiltonian Hr can be written in the standard 
way 

Hr = hur{a)a + l/2), (12) 

where 
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ih 

Pr = -Ti 



huJr 



4 ^ Ec,, 

Coupling term. The couphng term can also be written using the operators a, 



C,c 



A 



huJr 



E, 



C,r 



{2n - Ng) (a^ - a 



(13) 
(14) 



(15) 



Note that the characteristic coupling energy is Ec = yhujr/ EcrEccf^- 

A general analysis of the Hamiltonian H is beyond the scope of the present paper and 
will be presented elsewhere^. In the next section we will discuss an explicit matrix form 
of H, which can be obtained under certain simplifying conditions which are nevertheless 
experimentally relevant. 



IV. EIGENSTATES AND ENTANGLEMENT 

Throughout this section we will work in the zero-temperature limit. We are interested 
in the situation Ng ~ 1, such that we have to consider the charge qubit states |— ) and 
1+) only. Furthermore, as far as the resonator is concerned, we will consider huj. = Ej 
and work only with the ground state |0) and the first excited state In the limit of 
weak coupling, Ec <^ hur, the coupled system can be characterized by the four basis states 
|— , 0), | — , 1), |+, 0), |+, 1). The Hamiltonian matrix for this low-energy subspace reads 



H 



( Eo 


tEp 





-2ia(3Ec \ 


-tE(3 


El 


2iaf3Ec 








-2iaf3Ec 


E2 


iEa 


V 2iaPEc 





-iEa 


Es J 



(16) 



which is a hermitian matrix describing the two-level system coupled to the lowest states of the 
resonator. Here, Eq = E_ + hUr/2, Ei = E^ + 3huJr/2, E2 = E^ + hur/2, E^ = E^ + 3hur/2, 
Ea = Ec{2a^ - Ng), and Ep = Ec{2(3^ - Ng). 

Suppose that the system has been prepared in the state |'?/'(t = 0)) = |+,0) at time 
t = 0. This situation can be achieved by a suitable manipulation of the gate voltage Vg 
at times prior to t = oilii. At times t > 0, the time evolution of \ip{t)) describing the 
system is governed by the Hamiltonian (p^B]). We keep Vg fixed such that iV^ = 1 at t > 0. 
Thus we have = = 1/2, and hence E^ = Ep = 0. Moreover, as huUr = Ej, we have 
El = E2 = Ecj + Ej = E: without coupling, the state |+, 0) would be degenerate with the 
state | — , 1). Thus the Hamiltonian takes the simple form 



H 



f Eo -lEc \ 
E iEc 
-iEc E 
\tEc E3 J 



(17) 
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Note in particular that the state |+, 0) couples to the state| — , 1); as a result, the degeneracy 
between them is lifted and the states become entangled. The precise form of the entangle- 
ment is governed by the central 2x2 block of the matrix (p!7D. The eigenstates of this block 
are 

IXi) = [\-,l)+t\+,0)]/V2, 
\X2) = [\-,l)-^\+,0)]/V2, 

corresponding to the eigen energies E — and E + E^., respectively. These two excited 
eigenstates thus correspond to a maximum entanglement of charge qubit and resonator 
states, induced by the capacitive coupling between them. 
The time evolution of is given by 

\m) = [e-^^^-^^^^/'^lxi) - e-^(^-^^^)*/^|x2)] . (18) 
v22 '- 

We deduce that the state lipit)) oscillates coherently between |— ,1) and |+,0). In fact, 
these so-called quantum Rabi oscillations can be interpreted as the spontaneous emission 
and re-absorption of one excitation quantum by the resonator. An interesting quantity is 
the probability -Pi(t) to find the harmonic oscillator in the state |1) after a certain time t. 
This probability shows Rabi oscillations as a function of t with frequency 2Ec/h, 

P,{t) = |(l,-|^(t))|2 = i[l-cos(2i?,tA)]. (19) 

Since these Rabi oscillations are characteristic for the entanglement realized in the system, 
their measurement would provide direct evidence of the presence of the entangled states 
and 1x2)- We will discuss the feasibility of such a measurement in the next Section. 



V. DISCUSSION 

For the numerical estimates presented below we will consider parameters related to a 
typical aluminium superconducting circuitS'i. For the Josephson charge qubit we have cho- 
sen the following characteristics: Ej = 26.1/ieV, Ecj = 70/ieV, A = 240/ieV. As for the 
resonator, we take = 90pH and Ec\r = 12neV, as a result hur = 26.1/ieV. Finally, 
the coupling capacitance is chosen to be of the same order of Cj, Cc = 0.5fF, yielding 
Ec = 256neV. Note that the coupling energy is indeed much smaller than the Josephson 
energy, which in turn is equal to the excitation energy of the resonator. 

Using the above paramaters, we have plotted -Pi(t), Eq. ([T9|), as a function of time in 
Fig. H. We clearly see the Rabi oscillations with perdiodicity T^abi ~8 ns. 

In order to be able to observe these oscillations, we need to satisfy various conditions. 
First of all, in order to avoid the presence of quasiparticles, the temperature must be much 
lower than the gap A. Secondly, it is necessary to have a decoherence time which is longer 
than Tpiabi- In our system, the decoherence time will be the shorter of the lifetime of 
the excitated state of the resonator and the decoherence time Tqubit of the charge qubit. A 
Q-factor of about 500 is quite realistic for a superconducting LC resonator. This yields a 
lifetime > 10 ns. As for the qubit, the experiment by Nakamurai has indicated that 
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"Tqubit > 2ns. This lower bound is essentially the decoherence time associated with the cou- 
pling to the measuring probe. In principle, this time can be improved upon by modifications 
of the experimental set-up. Theoretical estimates! show that a time Tqubit ~ 100 ns should 
be feasable. Finally, a measurement of the number of excitations should be performed on the 
resonator. This can be done, e.g., along the lines of Ref.0, where the discrete, oscillator-like 
energy levels of an underdamped Josephson junction were measured. 
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FIG. 1. The quantum circuit. 
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FIG. 2. Probability Pi{t) to find the system, prepared in the state |+,0) at t = 0, in the state 
— , 1) after a time t. 
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